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Анализ и прогноз поведения судна 
в экстремальной ситуации 
на основе нечеткой системы знаний 


Обсуждается проблема повышения надежности принятия решений по обеспечению безопасности эксплуата- 
ции судна в экстремальной ситуации. Контроль и прогноз развития ситуации реализуется на основе нечеткой 
системы знаний. Моделирование осуществлено с использованием метода функционала действия. Особое 
внимание обращено на поведение судна при периодическом ухудшении остойчивости на волнении. 


Введение 


Исследование остойчивости судна на попутном волнении при периодическом 
снижении остойчивости — одно из важных направлений в общей проблеме обеспече- 
ния безопасности мореплавания [1-10]. Критические ситуации достаточно подробно 
изучены с помощью теоретических и экспериментальных методов. [6-8]. Наряду с 
ситуациями, которые стали типичными при изучении остойчивости на попутном вол- 
нении, существует еще одна ситуация, связанная с непрерывным ухудшением на- 
чальной остойчивости. Эффект снижения остойчивости в наибольшей степени прояв- 
ляется у судов с 5-образной диаграммой остойчивости. Наличие такой диаграммы 
при воздействии крупных пакетов нерегулярных волн вызывает резкое ухудшение 
начальной метацентрической высоты, что приводит к валкости судна. Известны слу- 
чаи тяжелых аварий, связанных с резким ухудшением остойчивости на волнении. 
Так, например, в работе [2] описана авария теплохода «Поронайск» на попутном вол- 
нении. Возникновению аварийной ситуации предшествовала сильная бортовая качка 
с амплитудой 30 — 35°, вызванная периодическим снижением начальной остойчивости 
(начальная метацентрическая высота уменьшалась до — 0,30 м и более при воздейст- 
вии крутых волн). 

В экспериментах с моделью грузового теплохода «Поронайск» на нерегулярном 
волнении [6] было установлено резкое ухудшение остойчивости, приводившее к вал- 
кости судна во время прохождения вершины волны через его мидель. Это способст- 
вовало возникновению бортовой качки большой амплитуды (до 30° и более), значи- 
тельно превышающей величину резонансной амплитуды (около 23°) при положении 
судна лагом к волнению. Из фрагмента записи изменения во времени начальной 
метацентрической высоты теплохода «Поронайск» при 8-балльном волнении (рис. 1) 
очевидно, насколько сильные колебания претерпевает величина в этих условиях. 
Следует отметить, что чрезмерный крен на попутном волнении отмечался и на дру- 
гих морских судах, для которых характерно резкое уменьшение начальной метацентри- 
ческой высоты, которое обусловило валкость при прохождении вершины волны через 
мидель судна [2], [6]. 


434 «Искусственный интеллект» 372009 


Анализ и прогноз поведения судна в экстремальной ситуации... ан 


( П И ПИПИТ 
05 1 


0 НИ | ИННА АМЕР 
| | | | 
=05 в И ЧИН тв Ц 1 УИ | | 
” 0 20 #0 60 80 АА 


Рисунок 1 — Изменение начальной метацентрической высоты при движении судна 
на нерегулярном попутном волнении 


В докладе обсуждаются особенности построения нечеткой системы знаний 
бортовой интеллектуальной системы (ИС), функционирующей на основе данных динами- 
ческих измерений, методов математического моделирования и структурированной базы 
знаний. При построении нечеткой системы знаний использованы результаты матема- 
тического моделирования динамики судна при периодическом снижении остойчивости 
на попутном волнении. Анализ проведен для судна с 3-образной диаграммой остойчивос- 
ти на основе дифференциального уравнения бортовой качки, включающего нелиней- 
ную пространственную функцию восстанавливающего момента на волнении. Основное 
внимание уделено анализу динамических картин поведения судна при различном 
уровне внешних возмущений. 


1. Особенности экстремальной ситуации 


Динамика судна на волнении описывается нелинейным дифференциальным 
уравнением [6]: 
(1х + цеоо)0"' + Мь(0’) + М(0,ф,5) = М» (0, (1) 


где (Зх + нев)0”, Мь(0’), М(60,ф5), М‚(® — функции, описывающие судна как динами- 
ческую систему (инерционные, демпфирующие, восстанавливающие и возмущающие 
компоненты). 

Наиболее сложной функцией в (1) является восстанавливающая компонента, 
входящая в дифференциальное уравнение бортовой качки. Восстанавливающая ком- 
понента отличается существенной нелинейностью, сложностью и многозначностью. Не- 
прерывно изменяясь во времени и пространстве, эта функция в значительной степени 
определяет результат интегрирования системы (1) при изучении физических картин 
взаимодействия судна с внешней средой. 

Математическое описание нелинейной пространственной функции восстанавлива- 
ющего момента на волнении представляется формулой [6]: 


М» =М(0,ф,0 =Б[(0,Ф) + АК0, 6) со$(о,1-=)]; (2) 
(0,ф) = 0,5 [(0,Ф)„ +1(0,Ф)ш |, АКО,ф) = 0,5 [4К0,Ф) м +АКО,Ф) ш |; 
Ми= Ф(0,фь)=р [(0,Ф, ий 


где А/9,ф)тах и АК@,Ф) ты — экстремальные значения приращений плеч остойчивости, 
соответствующие положению судна на подошве и вершине волны при различных 
курсовых углах ф; [(9,ф,{) — плечо восстанавливающего момента, определяемого для 
различных моментов времени, # — фаза: = =Ои 2л -— судно на подошве волны; & = д/2 — 
на переднем склоне; = = л — на вершине волны; & = 3/2 т — на заднем склоне. 
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Общее выражение функции приращения плеча остойчивости А/„(В»/А,0,Ф), опре- 
деляющее влияние волнения при различных параметрах формы корпуса и числа Фруда, 
в диапазоне углов крена до 60° включительно и ф = 0-:45° имеет вид: 


6 
Ар, (№, /1,0,Ф) = 814 "= бло, +547. 6.0) + 
т=1 (3) 


8 3 
+ В.В (0, Ф, ) +», СЕ, (0, Фу» 1, 
п=Т р= 


6 8 3 
Фр, /1,0,6,), >. А, (0,6, ), >. В,Е,(0,ф,), >.С,Е,(@,фл). 
т п=1 р= 

Входящие в модель (2) функции характеризуют данные стандартных моделей и 
суммы поправок на влияние линейных, квадратичных и кубических членов разложения 
ряда Тейлора на величину приращения плеча остойчивости, определяемого с учетом 
дифракционной и интерференционной составляющих при движении судна с про- 
извольным курсовым углом на волнении. Геометрическая интерпретация функции 
М(0,$,5) в виде сложной пространственной поверхности дана на рис. 2. 


М(0,</1) 


мМ(9,Ф,0 


к 
р. <» 
724% \ 


Рисунок 2 — Изменение нелинейной функции, описывающей восстанавливающий 
момент на волнении: сплошные кривые — мгновенные диаграммы остойчивости; 
пунктир — временные кривые (сечения поверхности при 9 = соп$0); 
штрих-пунктир — диаграмма остойчивости на тихой воде 


Вычислив приращения А/, в необходимом диапазоне углов крена, нетрудно 
установить плечи остойчивости /, и построить соответствующие диаграммы для 
положения судна на вершине и подошве волны по формуле /л = [+А/,, где / — плечо 
остойчивости, полученное обычным расчетом на тихой воде. При расчете на ЭВМ 
исходную формулу для /, удобно записать в виде 


И 
1, =1+ 8, +У Ая } (4) 
1=1 


где значения А; принимаются равными: 
А!* = ИВ - 4.82, А>* = В/Т - 2.67, Аз* = Н/Т - 1.3, 
А4* =к- 0.7, А5* =ф- 0.692, А‹* = Ег- 0.28, (5) 
Аз=В1=Аг, Аз=В›= А», Аэ= В: = Аз”, Аю=Ву= Аз, А =В5= Аб, А12= Вв= А! А», 
Аз = Ву= ДА, А14= Вз= А, Де, Аб = С1= А», Ав= С›= Аз’, Ал7= Сз= Аз; 
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где 1, |, К- счетчики циклов: 1 определяет номер поправки на отклонение характеристик 
формы корпуса и числа Фруда от данных стандартной модели, 0 - углы крена, 
ф- курсовые углы. 


2. Нечеткая формальная система знаний 


При функционировании ИС обеспечивается контроль ситуации и выработка 
практических рекомендаций для принятия решений по выбору оптимального режима 
движения в зависимости от интенсивности внешних возмущений. Процесс управления 
организуется в рамках концепции «мягких вычислений» [10]. Управляющий алгоритм 
представляется набором логических правил «Е — ТНЕМ» [1], [8]. В качестве текущей 
информации используются показания датчиков колебательного движения судна при 
бортовой и килевой качке. 

В бортовых ИС находят применение модели нечеткого вывода, основанные на 
различных формальных конструкциях, обеспечивающих разрешение нечетких логи- 
ческих правил. При реализации измерительной и экспертной информации в ИС при- 
меняются правила нечеткого вывода, использующие принцип обобщения Л. Заде для 
нечеткого вывода, в котором высказывания моделируются нечеткими множествами, 
позволяющими вывести нечеткое заключение. Основная проблема реализации нечеткого 
управления в бортовых ИС заключается в построении функций принадлежности (ФП), 
особенно в условной части логического правила, содержащего несколько нечетких 
переменных. Построение ФП осуществляется на основе обработки измерительной и 
экспертной информации, а также материалов конкретного проекта судна. При этом 
реализуются два подхода. Первый подход состоит в формализации нечетких знаний 
в процессе разработки ИС. Он основан на традиционном представлении экспертных 
знаний и предусматривает построение ФП базовых правил, «зашитых» в программный 
комплекс. Второй подход реализуется в процессе функционирования ИС и связан с 
построением текущих значений ФП на основе данных динамических измерений. 
Алгоритмы и программное обеспечение построения ФП основаны на использовании 
статистических методов и нейросетевых технологий [7]. 

Композиционная модель нечеткого вывода описывает связь всех возможных 
состояний логической системы с управляющими воздействиями и формально запи- 
сывается в виде (Х,В,У), где Х = {хь,...,Хи}; У = ф$уь....,Уш} — базовые множества, на 
которых заданы входы А; и выходы В; нечеткой системы; К. — нечеткое соответствие 
«вход — выход», которое строится на основе формализации нечетких стратегий. 
Принятие решений при функционировании нечеткой системы знаний осуществляется 
на основе модели [1], [8]: 


«ША, — Феп Вт, ее ЁА»>› Феп В....., е1зе № Ах Феп Вх». (6) 


Здесь Ау1,...,Ах — нечеткие подмножества, определенные на базовом множестве 
Х,аВ.,....,Вм — нечеткие подмножества из базового множества У. 

Антецедент импликации Х представляет собой проверку выполнения условий, 
<Х, в А, & Х) 18 А>, сформулированных на основе соотношений (11). Консеквент 
импликации У(7,) определяет возникновение экстремальной ситуации и выработку 
соответствующих практических рекомендаций, в том числе и выполнение операций 
по моделированию и визуализации ситуации. Функции принадлежности, характеризую- 
щие условия возникновения экстремальной ситуации, представлены на рис. 3. Здесь 
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значения ©! и 6› характеризуют диапазон перехода к полной потере начальной остой- 
чивости, АБ! и ЛЁ› — соответствующие значения приращений метацентрической высоты 
в момент прохождения вершины волны через мидель судна. 


Рисунок 3 — Функции принадлежности антецедентной части логического правила, 
описывающего возникновение рассматриваемой экстремальной ситуации 


Для решения задач принятия решений на основе априорной информации 
используется нечеткая формальная система, обеспечивающая принятие решений при 
выборе альтернатив [8] в рамках концепции мягких вычислений ой сотрийпе [10]. 
Формирование такой системы осуществлено на основе принципов обработки инфор- 
мации в мультипроцессорной вычислительной среде [7]. 

Нечеткая формальная система разработана с учетом динамики изменения пред- 
метной области в процессе функционирования ИС. Структура расширенной нечеткой 
системы (Ежепаеа Еитху Зу\бет — ЕЕ) представляется в виде матрицы переходов, 
строками которой являются пары [3]: 


(ЕЕЗ,, «Го 1= Е] и И (7) 


где ЕЕЗ1 — расширенная нечеткая формальная система; <Т:,$}> — пара, задающая мак- 
симинную операцию, используемую в основе композиционного правила вывода (0боб- 
щенный тоЧиб ропепз); Т‚— оператор минимума (Т-норма); 5] — оператор максимума 
(Т-конорма). 

Пересечение строки ЕЕ$к и столбца ЕЁЕ$, матрицы переходов обозначается 
символом «Х» и указывает на возможность выполнения максиминной операции, 
определяемой парой <Т,,5}>. Эта операция сменяется на максиминную операцию, за- 
данную парой <Ть,5„>. Условия переходов определяются заранее при разработке не- 
четкой формальной системы. 

При разработке нечеткой модели знаний учитывалась возможность переходов 
от одной расширенной нечеткой формальной системы к нескольким подобным сис- 
темам, которые могут создаваться и удаляться в процессе эксплуатации на основе 
моделирования динамики судна в рассматриваемой экстремальной ситуации. В этом 
случае двумерная матрица переходов заменяется на трехмерный куб переходов, 
располагаясь вдоль осей Х,У,/, соответственно: 

— Х- множества {ЕЕ$1,...,НЕЗ и}; 
— У- множество Т-норм; 
— Д- множество Т-конорм. 

При функционировании ИС происходит накопление новой информации и 
корректировка базы знаний. Динамическая база знаний совершенствуется и расширяется 
за счет включения новых моделей. 
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3. Математическое моделирование динамики судна 
в условиях непрерывного изменения поперечной 
метацентрической высоты судна на волнении 


В процессе моделирования исследовалось влияние параметров волнения на 
изменение мгновенных ординат диаграммы остойчивости /[(0,1) в условиях развития 
шторма. Экстремальные значения ординат этих диаграмм устанавливались для различ- 
ных значений крутизны волны (/А/Вху), величина которой постепенно увеличивалась 


по мере прохождения пакетов волн в моменты времени ц, Ц...., &. Приращение 
начальной метацентрической высоты ДВ вычислялось по формуле: 
ДВ = Бо г (АБ) инь (8) 


где (АБ) — мгновенное значение метацентрической высоты при положении судна 
на вершине волны, определенное на основе математической модели (2) - (4). 


Рисунок 4 — Кривые изменения метацентрической высоты для различного 
уровня внешних возмущений: А — случай нормальной эксплуатации 
(явление валкости не проявляется); В — типичная ситуация, когда при 
воздействии пакета волн появляются отдельные отрицательные значения 
метацентрической высоты; С — предельный случай, соответствующий 
полному ухудшению начальной остойчивости 


Сценарии исследуемых ситуаций предусматривали также различные случаи 
поведения судна при смещении груза вследствие резкого снижения начальной 
остойчивости. Особое внимание обращалось на трансформацию диаграммы остойчи- 
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вости при прохождении вершины волны через мидель судна в условиях интенсив- 
ных волновых возмущений. На основе результатов моделирования сформулиро- 
ван критериальный базис, обеспечивающий предотвращение возникновения опас- 
ных ситуаций вследствие резкого уменьшения начальной метацентрической высоты 
при движении судна на попутном волнении. Критерии остойчивости использованы 
при разработке базы знаний бортовой ИС контроля динамики судна на волнении с 
учетом расширения нечеткой формальной системы. 

Результаты моделирования экстремальных ситуаций при движении судна на нере- 
гулярном попутном волнении представлены на рис. 4. Здесь отобраны три харак- 
терные ситуации, которые могут наблюдаться в условиях эксплуатации исследуемого 
судна. Они представляют собой нестационарные режимы колебаний метацентричес- 
кой высоты под воздействием случайных возмущений (прохождение пакетов волн с 
непрерывно увеличивающейся крутизной). Штриховкой на этих рисунках указана 
область отрицательных значений метацентрической высоты. 

Результаты моделирования показали, что возникновение «раскачивающего 
эффекта», обусловленного валкостью судна, обнаруживается в ситуациях, когда 
величина начальной метацентрической высоты судна с З-образной диаграммой 
остойчивости и соотношения параметров развивающегося волнения, представленного в 
виде пакетов нерегулярных волн, достигают следующих значений: 


№ < 0,3 м, ^/ Г = 0,8 + 1,3, Ву /А, > 0,08. (8) 


Оценка устойчивости колебательного движения судна в условиях периодического 
ухудшения начальной остойчивости осуществлялась на основе метода функционала 
действия [10]. Разработанный алгоритм позволяет вычислить вероятность опрокидывания 
по мере развития колебаний судна, вызванных появлением отрицательных значений 
метацентрической высоты. 

Стохастический алгоритм, реализующий последовательность операций при 
выполнении процедур метода функционала действия, представлен на рис. 5. Здесь 
Мх0”), Мь(0’), М(0,Ф.0, М(@® — инерционно-демпфирующие, восстанавливающие и 
возмущающие компоненты математической модели; & — компонента гауссовского 
случайного процесса (®) >В”, формирующий фильтр которого определяется параметрами 
доминирующей частоты © и шириной спектра возмущений Е; У(®) — «белый шум» 
малой интенсивности; \ — винеровский случайный процесс; = — малый параметр; 9О — 
область неустойчивости. 

Анализ устойчивости колебательных режимов качки корабля на волнении в 
условиях периодического снижения начальной метацентрической высоты позволил 
установить ряд характерных закономерностей динамики исследуемой нелинейной 
системы. Эти закономерности определяются особенностями групповой структуры набе- 
гающих волн и интенсивностью колебательных режимов. Результаты вычислитель- 
ных экспериментов позволяют построить и проанализировать траектории динамичес- 
кой системы на фазовой плоскости [6], [7], [9]. Наиболее типичными колебательными 
режимами являются случаи возникновения аттракторов, описывающих динамику 
одностабильной, бистабильной и трехстабильной систем [8], а также характерные 
временные кривые и одномерные отображения из точек, соответствующих максималь- 
ным значениям амплитуды для (п + 1) и п-го интервалов времени для аттракторных 
множеств [4], [5]. 
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Анализ и прогноз поведения судна в экстремальной ситуации... ан 


Исходная модель 
(1х +нов)0”' + Мь(60’) + М(0,0,5) = М@ 
Формальная модель 
У) =Ы((0))) + 6(5(0))(@); У) = (0",0',6,,6,)' ® В" 
Уд =о’, =>0; 61) = Е (+ Су, 6(0) = 6, 


==] 1 6- 1=Ьа&, =. (0? +Е?) -2Е; 
| 8. | Х=-2Е, Г, =-6* -Е?. 

® Реализация алгоритма 
\У+=[ю,Т] найти & и управление ПО из заданного класса непрерывных 
по [ю,Г] функций, обеспечить граничные условия и У(0) = ус, У()ЕДО 
и минимизировать функционал 


1 Т 
по =5 [у Фа 
2 0 
® Вероятность опрокидывания 
Р(у(еБО) определяется асимптотическим соотношением 


та” ш Р{у(Р) Е ОО} =-шщ 7 (1) 


» Графическая интерпретация 
Результаты моделирования - в виде зависимостей 0 (1), Р(+) 
в области устойчивости © = Ё(у,ф) 


Рисунок 5 — Алгоритм, реализующий последовательность операций 
на основе метода функционала действия 


Интересная закономерность, установленная в процессе вычислительного экспери- 
мента, заключается в моделировании пространственно-временного поведения судна 
на базе центральной теоремы синергетики, определяющей динамику сложной систе- 
мы как взаимосвязь растущих и затухающих ее конфигураций с параметрами поряд- 
ка [4], [5]. Установлено, что все исследуемые пространственно-временные состояния 
системы управляются параметрами порядка (принцип подчинения). Переход к пара- 
метрам порядка в сложной системе приводит к сжатию информации и к значитель- 
ному сокращению числа степеней свободы. При этом параметры порядка определяют 
движение элементов системы, которые, в свою очередь, совместно (кооперативно) 
характеризуют действие параметров порядка (принцип круговой причинности) [4]. 


Заключение 


Таким образом, технология обработки информации в бортовых ИС при анализе 
и прогнозе динамики судна в условиях непрерывного изменения метацентрической 
высоты имеет следующие преимущества: 

1. Расширение традиционных подходов к обработке информации в сложных 
динамических средах, дополнение их новыми методами, моделями и алгоритмами 
поддержки принятия решений по управлению судном в сложной гидрометеорологи- 
ческой обстановке. 

2. Учет неопределенности и неполноты исходной информации при контроле 
ситуации и выборе решений по управлению судном в условиях непрерывного изменения 
метацентрической высоты. 
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3. Развитие внутреннего потенциала теории принятия решений на основе принципа 
конкуренции, позволяющего на базе анализа альтернатив выбирать предпочтительную 
вычислительную технологию с использованием методов классической математики и 
нечетких моделей. 
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Ве 


Ю.Л. Нечаев, Хейн Тун 

Анал 1 прогноз поведнки судна в екстремальнй ситуаций на основ! нечтко: системи знань 
Обговорюеться проблема шдвищення надйност! прийняття рипень 13 забезпечення безпеки експлуатаци 
судна в екстремальний ситуаций. Контроль 1 прогноз розвитку ситуацй реалзуеться на основ! нечитко! системи 
знань. Моделювання здйснено з використанням методу функцюнала дй. Особливу увагу звернено на 
повед1нку судна при перодичному погршенн! остиност! на хвилюванн!. 


Уи.1. Месйаеу, Нет Тип 

'Тве Апа]у$5 ап@ Еогеса$ о ЗШр Вепау!ог ш Ехбете ЗИиайоп оп Че Ва$5 оЁ Еи7лу Кпо\едое Зу$ет 
ТБе ргоет оЁ 1шсгеазе оЁ тепабИИу о ассер{апсе оЁ Фе 4ес1510п$ оп а за у оЁ зв1р орегайоп ш ап 
ехгете зйаайоп 1$ 915си5зе4. ТВе сопёго| ап Ююгеса${ оЁ 4еуеюртеп{ оЁа зИлайоп 15 геаП7е4 оп е Ъаз1$ 
ОР 7ту Кпо\Ле4дее зузбет. ТВе тодепп» 15$ саглей ой у зе оРа Фе ВпсНопа| асНоп тефо4 ТБе 
зресла! айепНоп 1$ шуеце4 оп $51р Беваутог а{ репо с еепогайоп забИиу оп \ауез. 
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